Nonadiabatic conditional geometric phase shift with NMR 
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Abstract 

A conditional geometric phase shift gate, which is fault tolerant to certain types 
of errors due to its geometric nature, was realized recently via nuclear magnetic 
resonance(NMR) under adiabatic conditions. By the adiabatic requirement, the 
result is inexact unless the Hamiltonian changes extremely slowly. However, in 
quantum computation, everything must be completed within the decoherence 
time. High running speed of every gate in quantum computation is demanded 
because the power of a quantum computer can be exponentially proportional 
to the maximum number of logic gate operations that can be made sequentially 
within the decoherence time. The adiabatic condition makes any fast conditional 
Berry phase(cyclic adiabatic geometric phase) shift gate impossible. Here we 
show that by using a newly designed sequence of simple operations with an 
additional vertical magnetic field, the conditional geometric phase shift gate can 
be run nonadiabatically. Therefore geometric quantum computation can be done 
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at the same rate as usual quantum computation. 

A fault-tolerant quantum logic gate jlj is the central issue in realizing the basic constituents 
of a quantum information processor. Quantum computation via the controlled geometric phase 
shift provides a nice scenario to this purpose. Due to its geometric property, a geometric 
phase |,1 shift can be robust with respect to certain types of operational errors. In particular, 
suppose that the spin(qubit) undergoes a random fluctuation about its path in the evolution. 
Then the final value of the geometric phase shift will not be affected provided that the random 
fluctuation does not change the total area ||. 

Recently, it was reported @J^-3 that the conditional Berry phase( adiabatic cyclic geometric 
phase) shift can be used in quantum computation. In particular, in ref |J, an experiment was 
done with NMR [|10| -|14|| under the adiabatic condition. 

Demanding on both running speed and precision of every gate in a quantum computer is 
quite high. It has been reported recently that ||, due to the requirement of adiabatic condition, 
faster running speed causes severe distortions to the results [f2| . The adiabatic condition makes 
the fast speed and high precision conflict each other in the conditional Berry phase shift gate. 
Increasing the running speed of the logic gate can exponentially increase the power of a quantum 
computation such as quantum factorization [|15|]. It seems that the geometric phase shift gate 
will only be practical if it can run at a speed comparable to that of the usual quantum gate. 
Therefore one is tempted to find an easy way to make the conditional geometrical phase shift 
through nonadiabatic state evolution. In this letter, we present an easy scheme to make the 
geometric phase shift nonadiabatically with NMR. We start from the exact control of the state 
evolution on a cone. 

Exact state evolution on the cone. It is well known that a spin half nucleus can gain 
a geometric phase shift in the conical evolution. We now assume that the initial spin state 
is on the cone and we will demonstrate a method to make it evolve nonadiabatically. The 
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Hamiltonian for a spin in a constantly rotating magnetic field [] is 

H(t) = [u (T z + U 1 <T x (t)]/2. (1) 

( o e'^ 1 ^ 



Here a x (t) 



e iit Q 



ujq is the amplitude of vertical field and U\ is the amplitude of 



the horizontal field, which is rotating around z— axis in the constant angular speed 7. The 
initial state \ipo) is an eigenstate of H = H(0). Explicitly, |^o) — cos|| |) + sin|| |), 



cos 9 = ujq/Juq + u)\ . To know the cause of the state distortion, we study the time evolution 
first. (The similar idea has been used for the adiabatic rotational splitting with NQR previously 
|]IE|| .) The time evolution operator U generated by the Hamiltonian H is determined by the 
time-dependent Schrodinger equation 

i^U = H(t)U . (2) 

Solving this equation we obtain 

U (t) = e -i"f^ e -i(Ho-l^/2)t_ (3) 

In particular, at time t = t = 2-71-/7, when the external field completes a 2tt rotation, the state 
is evolved to 

tP(r) =e~ ln - i{H °-^ /2)T \^ }- (4) 

The adiabatic approximation assumes that at time r, the state completes a cyclic evolution 
provided that 7 is small. However, due to the decoherence time limitation, 7 cannot be too 
small. The non-zero 7 distorts the state in the evolution and makes it noncyclic at time r . This 
non-zero 7 in the adiabatic approximation causes two types of errors. One is that the noncyclic 



We shall extend it to the case of arbitrary time-dependent rotating speed in the end of this section. 
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completion of the evolution at time r will cause further errors in the succeeding operations in 
the sequence. (A sequence of operations is used to remove the dynamic phase in realizing the 
geometric quantum gate ||.) The other is that the geometric phase acquired over period r is 
not — 7r(l — cos 6), as for the ideal adiabatic cyclic evolution. 

Now we give an easy way to exactly control the state evolution on the cone. Here the 
external field can be rotated around z— axis arbitrarily fast. We use \ipo), the eigenstate of H 
for the initial state. We switch on a static vertical magnetic field u z = 7 while the external 
field is rotated around z— axis. In this way, the new time-dependent Hamiltonian is H\y{t) = 
I [{ujq + "i)a z + uia x (t)]. The time evolution operator Uw(t) generated by this Hamiltonian 
satisfies the Schrodinger equation 

i^U w (t) = H w (t)U w (t) (5) 

with the boundary condition Uw{0) — 1- The state at time t is related to the state at time 
by \ijj(t)) = Uw(t)\ipQ), Vo — ^(O). Denoting R = e n * CTz//2 we obtain the following equivalent 
equation 

td(RU w )/dt = Ho(RU w ). (6) 

Noting that H is time-independent therefore we have RUw{t) = e~ lHot . This is equivalent to 

U w (t) = e -^ ta ^ 2 e- iHot . (7) 

Consequently, at any time t, state \ip(t)) = e~ lXt e~ l/ytaz ^ 2 \ipo) is exactly the instantaneous eigen- 
state of Hamiltonian H(t), where A is the eigenvalue of Hq for eigenstate |^o)- In particular, 
at time r, = e~ m ~ lXT \il)Q) . It only differs to j^o) by a phase factor. Note that the time 

evolution operator Uw(t) here is generated by the Hamiltonian Hw(t) instead of H(t). 

Here the additional vertical field 7 plays an important role. Without this field, the time 
evolution operator generated by H(t) is 
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U(t) = e -nta z /2 e -iH lt (g) 

and Hi = H — / ~fa z /2. Obviously, this time evolution operator U(t) will distort a qubit with 
the initial state \ipo) in the evolution. But this U(t) can exactly control the qubit with the 
initial state the eigenstate of H±. If initially we set the spin state to l^i), then the spin 
will evolve exactly on its cone without the additional field 7. So, we have two ways to control 
the spin evolution exactly. We can use the additional field 7, if the initial state is set to be 
l^o)- Alternatively, we can also set the initial state to be and then we need not add any 
additional field when the external field is rotated. In this letter, we adopt the former one, i.e. 
we set the initial spin state to |^o). Q 

We have assumed above that the field rotates at a constant speed. Actually, we can easily 
modify the above scheme for arbitrary time-dependent rotating speed 7 a (t), Jq 7 a (t)<it = 2n. 
In this case, we need only change the term jt in H w (t), Uw(t), R and ip{t) into JaJait^dt' 
accordingly. Consequently, the the additional vertical field is now a time-dependent field u z (t) = 
7 a (t) instead of a static field. This extension is important in case it is difficult to rotate the 
field (or the fictitious field 0) in a constant angular speed. Punctual results can be obtained 
here through the exact feedback system where the value of additional vertical field is always 
instantaneously equal to the angular velocity of the rotating field. 

Thus we see, by adding an additional magnetic field that is equal to the rotating frequency 
of the external field, we do get the exact result. Obviously if we rotate the field inversely, 
additional vertical field should be in the inverse direction(— z) accordingly. 

NMR system and the rotational framework. Consider the interacting nucleus spin 

2 Our main motivation is to give the nonadiabatic quantum gate in the most general case, i.e., we 
shall allow the time-dependent rotation of horizontal field uj\. For this purpose, only the method using 
additional vertical field will work. 
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pair(spin a and spin b) in the NMR quantum computation [p~0| [T3|1 . If there is no horizontal 



field the Hamiltonian for the two qubit system is Hi = \{u a a za + 0Jbcr z b + Jo za ■ <r z b), where 
uj a (ujb) is the resonance frequency for spin a(b) in a very strong static magnetic field(e.g. uj a can 

(l 

be 500MHz , J is the interacting constant between nuclei and a za = a z b = a z = 

I - 1 / 

After adding a circularly polarized RF field in horizontal plane, the Hamiltonian for spin a in 
static framework is 



rr/ 1 1 / 1 

H = -lu (J z + -cu a a z + -u x 



e 
y e iJ - 1 j 



(9) 



Here u' a and u)\ are the angular frequency and amplitude of the RF field respectively, ojq = 
u a — uj' a ± J and the " ± " sign in front of J is dependent on the specific state of spin b, up 
or down respectively. To selectively manipulate spin a, we shall use the rotational framework 
that is rotating around z— axis in speed uj' a . We assume u' a close to u) a but obviously different 
from uJb- The Hamiltonian for spin a in the rotational framework is 

H a = R'H'R'- 1 + i(dR' /dt)R'- 1 = H (10) 

and R' = e tuJ ' a(7zt ^ 2 . Note here H is dependent on state of spin b through loq. In the rotational 
framework, if the horizontal field is rotated in the angular speed 7, the Hamiltonian for spin a 
is just Hit), as defined in eq(|l]). In the NMR system, we require \ip ) to be the eigenstate of H 
in rotational the framework. Previously Q, state \tpo) for spin a was produced adiabatically. 
In the following we demonstrate how to nonadiabatically produce the state \ipo). 
Creating the conditional initial state with NMR. In rotating the state around z— axis, 
we have started from state |-?/>o)> which is the eigenstate of H . State |^ ) must be created from 
the spin state | f) or | j). We need a nonadiabatic way to create state |^o) f° r s P m a - Note 
that we are working in the rotational framework. We denote 5 = u a — u' a . We can use the 
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following sequence of operations to create the conditional angle 8. (For simplicity we will call 
the following sequence as S operation later on.) We have S = 



J'0±(t c )) -> [S-t c 



<P? (11) 



Here all terms inside the [• • •] represent the Bloch sphere rotation angles caused by RF pulses. 
The superscripts indicate the axis the Bloch sphere is rotated around. J'(<p±(t c )) is the time 
evolution over period t c by the Hamiltonian |(5 ± J)a z . This evolution rotates spin a around 
z— axis for an angle ip± = (S ± J)t c . (Or (p± + tc, if spin a is down initially. For clarity we omit 
this case and alway assume spin a is initially up.) After this S operation, the angle between 
spin a and the z-axis is 9± = f — (</?' + v?±)(see Fig. 1). To ensure the state to be eigenstate of 
Hamiltonian H$ after S operation we require tan(<// + Jt c ) = and tan(</?' — Jt c /2) = 
simultaneously. This is equivalent to 

Jt c = (arctan ^ - arctan 12 

v wi uii ii (12) 

(ff = (arctan ^ + arctan ^) / 2 

From this constraint, given specific values of 5, J and U\, we can easily obtain the scaled time 
J ■ t c (Fig. 2) and the angle ip'(Fig. 3) in controlling the S operation. In particular, for 
8/ J = 1.058, the value adopted in the recent NMR experiment, the relation curves for Jt c vs 
uj\ and if' vs uj\ are shown in Fig. 4. 

The nonadiabatic conditional geometric phase shift. After the S operation, the state 
of spin a is an eigenstate of Hamiltonian H no matter whether spin b is up or down. After 
this S operation we can use the additional magnetic field to nonadiabatically control the state 
evolution on the cone. However, after the state completes a loop, the total phase shift includes 
both dynamic and geometric contribution [pjj. It is possible for the state to evolve on the 
dynamic phase path if we choose certain specific value of the rotating speed 7( and also the 
additional field 7). For this task, we require 



The negative sign in the right hand side indicates that tha additional field 7 is anti-parallel 
with the field u>o- By this setting, the total magnetic field in Hi is always "perpendicular" 
to the state vector expressed in the Block ball. One can easily show that the instantaneous 
dynamical phase 

< ^\U ] w {t)H{t)U w {t)\^ >= 0. (14) 

In case that 7 is time dependent, to remove the dynamical phase, we can use the time- 
dependent magnetic field Ui(t) and u> (t). We require 

= 

and 

w (f)M(f) = w o (0)/wi(0) = tan# (16) 

The two qubit case is a bit different from the single qubit case. Qubit b could be either up 
or down. We need choose the appropriate 7 value so that the instantaneous dynamic phase for 



qubit a is always zero no matter qubit b is up or down. For this purpose, besides equation QUID, 
we require 



w x = V5 2 - J 2 . (17) 

With this setting we propose the following scheme to make the nonadiabatic conditional geo- 
metric phase shift. 



S 



f uj z ' 



S~ l . (18) 



(., \ 



The term 



to 



represents doing operation C with an additional vertical field 7. C represents 



rotating the external field around z— axis for 2n in a uniform speed 7. We should choose 
the rotating direction so that additional field 7 is anti-parallel the field 8 ± J, (5 > J). The 
dynamical phase shift after the operations is 0. So the scheme raised here can remove the 
dynamic phase and retain only the geometric phase. It can be shown that the geometric phase 
acquired after the total sequence of operations is T + , r_, — r_ and — r + respectively for the 



four different initial state (| ft), | ff), | ff), | ff)), T± = — ir =|= 2n cos9±, cos9± 



8±J 



s±j 

26 ■ 



Concluding remark. To make the conditional geometric phase shift gate we need be able 
to control two elementary operations. One is to exactly control the cyclic state evolution on a 
cone; the other is to produce the initial state on the cone with an angle 6± conditional on the 
other bit. Previously 0, these two elementary operations were done adiabatically. We have 
shown that both of these two tasks can be done nonadiabatically. We have designed a new 
scheme for the nonadiabatic conditional geometric phase shift gate. This makes it possible to 
run the geometric quattum gate in a speed comparable to that of the normal quantum gate. 
The idea on nonadiabatic geometric phase shift gate demonstrated by the NMR system here 
should in principle also work for the other two level systems, such as Josephson Junction system 
H and the harmonic oscillator system [p,|IB], where the decoherence time can be much shorter. 
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Note added: This manuscript is a revised version of Ref. [19] with corrections made according 



to Ref. [ 20 1 . In this version, we removed all unpublished contents in the old versions therefore 
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the new version includes only the contents published in |T9| , |20| . For an extended version with 
all details, one may refer to |quant-ph / 1 08 1 1"T| pTf , which was posted to the arXiv one year ago. 



The original scheme appeared in Ref. (19| does not filter the dynamical phase. This was first 
pointed out by Tim Piller et al in August 2001. We then replaced the old multiloop scheme by 
a new single loop scheme with dynamical phase free path. This correction was first presented 
as the version 4 of this manuscript and also as |quant-ph/01081iTl in August 2001, and finalized 



in 2002 20 
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FIGURES 
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FIG. 1. Using 5 operation to produce the eigenstate of Hq nonadiabatically. Picture 

a shows the possible initial state for spin a. The up(down) arrow on the solid line represents the 
up(down) state for spin a while spin b is up. The up(down) arrow on the dashed line represents the 
up(down) state for spin a while spin b is down. In picture f, angle GO A is 9 + = arctan j^j and angle 
GOB is 6L = arctan f^. 
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FIG. 2. The time control in S operation. Parameter uj\/J varies from to 10, 5/ J varies 
from to 5. Vertical axis is for the scaled time Jt c . 




FIG. 3. The rotating angle control in S operation. Parameter lo±/J varies from to 10, 5/ J 
varies from to 5. Vertical axis is for the angle ip' . 
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FIG. 4. The time and rotating angle control in S operation for specific 5/ J value. 

Horizontal axis represents uj\/J. Vertical axis is for the scaled control Jt c (the solid line) or the angle 
(//(dashed line). Here 5/ J = 1.058, as for the experimental condition[2]. By this figure, given the 
specific values of wi/J, we can always find the corresponding point in the two curves thus we can take 
the suitable time control(for t c ) and rotation control(for ip') in the S operation. 
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